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Abstract 

In quantum information theory, von Neumann entropy plays an important role. The entropies can be obtained 
analytically only for a few states. In continuous variable system, even evaluating entropy numerically is not an easy 
task since the dimension is infinite. We develop the perturbation theory systematically for calculating von Neumann 
entropy of non-degenerate systems as well as degenerate systems. The result turns out to be a practical way of the 
expansion calculation of von Neumann entropy. 

1 Introduction 

In quantum information theory, von Neumann entropy of a state appears in many basic theorems such as quantum source 
coding 1] , quantum channel coding of classical information [5] [3] , quantum channel coding of quantum information [3] 
(SJ- The later two are problems of channel capacities, which are maxima of Holevo quantity and coherent information, 
respectively (usually the regulation procedure should be taken). The problems can be reduced to the calculation of von 
Neumann entropy. The basic method to obtain the von Neumann entropy of a state is to calculate its spectrum. The 
spectrum can seldom be obtained for a given state. Meanwhile numeric calculation of the spectrum may encounter 
problems for continuous variable system which has an infinite dimensional Hilbert space. Hence a perturbation theory 
for the von Neumann entropy of a state is needed. While quantum mechanics provides us the theory of lower order 
perturbation to the energy levels, here we will develop a systematical theory of the entropy perturbation up to any 
order of precision. 

2 The perturbation to the entropy of a non-degenerate system 
2.1 The perturbation of entropy via eigenvalues 

A non-degenerate quantum system is a density matrix po with its eigenvalues E n ^ E rn for all n ^ m, where n (or m) 
is the number (or vector number) specifying the quantum eigenstate. The usual quantum perturbation theory for the 
non-degenerate systems gives the first and second order eigenvalue perturbations 

E n X) = H nn, (1) 

tip = ( 2 ) 

when the perturbation to the density matrix is the matrix H, where H nm = (ip n \ H \ip m ) , with \ip n ) the eigenvector of 
the density operator po corresponding to the eigenvalue E n . As usual, the perturbed density matrix can be written 
as p = p + eH. Since p and po are density matrices, we have W — H,and TrH = 0. Thus we have ^ n En = 0,and 
it is clear that J2 n = 0- ^ ne total eigenvalue up to the second order is E n ^ = E n + eE$p + e 2 En K The entropy 
of the state up to the second order perturbation is S(p) — —Trplogp — ~^2 n En' > log E^ + o(e 2 ). Thus we have 
S(p) = S(p ) + + ie 2 %^ + o(e 2 ) with 



dS(po) 



ds 



]r# im log£„ (3) 



de 2 ^ E n 
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2.2 The expansion formula of the entropy 

For a > 0, we have loga = / °° ^^-j^pr- Similarly, for a positive operator A, we have log^4 = J °° A £~^ . [5] Thus 
log A-logB = - liniM^oo J M (^j- -w+i)dt for positive operators A and B. ds ^ o) = - lim e ^ \Tr[p\ogp-p logp ] = 
-lim E ^oiTr[p (logp-logpo)] -Tr[Hlog Po \, Note that ^ = £ E£l (-l) fc (4)*, we have 



dS(po) 



/ Tr{po(po + «)- 1 ^(po + t)- 1 }^-T'r[lflogpo]- (5) 
Jo 



Evaluating in the eigenbasis of po, it reads 

= E + t)~ 2H nndt - H nn ^gE n 

n n 

= -^2H nn \ogE„, (6) 

n 

where TrH — has been used. It coincides with Eq.([3]), the result derived from the perturbation of the eigenvalues. 

The second derivative of the entropy at po is d ^li ^ — — lim e ^o j?Tr [p + log p + + p_ log p_ — 2po log po] , where p± = 
po±eH. The derivative can be rewritten as Pi + P2, the two parts are Pi = — lim e _»o -^Trpo(logp + + logp-~2\ogpo) — 
2 J^Trpoipo+t^H (p +t)- 1 H(p +t)- 1 dt a ndP 2 =-lim £ ^ ±Tr#(logp + -lo g( L) - -2 TrH(p +t)- 1 H(p + 
t)~ 1 dt. With the eigenbasis of po, we have 

d 2 S( Po ) f M t|g„ m | 2 

d £ 2 £-7 (£ n + ^(£; m +t) 

\ ~« I H nn \ \ - \H nm \ E m E m 

= -E%^- 2 E #5r l0 ^ (7) 

Note that Eq.(|?J) and Eq. JTJl are strictly the same, and the non-degenerate condition is used. 

The n th derivative of the entropy can also be carried out. From the definition of the derivative, we have d = 
- hm^o i ELo S(po+ls) Q = P 3 +P 4 , with P 3 = - lim^o £ ^=0 Tr <°° lo g(Po+^) (£) (-1)' - n! / °° Tr Po ( P o+ 
ty^Hipo+t)- 1 }" dt for n > 0. P 4 = - lim^o pr ELo Tr[ei.Enog(po + fe)] Q (-1)' = -n! / °° Tr[H(p + i)" 1 ]" <ft 
for n > 1 ( see Appendix A). Thus we have 

' / "' V/,nl --(-l) n n! /"rr^Cpo + tJ-M^ + t)- 1 ]"}*. (8) 



efe" 

The Taylor expansion formula for the entropy of a state p = p$ + ei? will be 

" e n d n S( Po ) 



s(p) - E 



n ! tfe n 

n=0 



= 5(p ) -eTr[P log po] - E(- £ )" / Tr{t(p Q + ty l [H(p a + ty l ] n }dt. (9) 

n=2 ^0 

We calculate the third and the fourth order perturbation for the case of H nn = 0, the entropy up to fourth order 
perturbation is 

s( P ) = s(p ) - e 2 1 H tI lQ g En + g3 Q 3 - g4 (Q 41 + Q 42 + Q«) + °( g4 ) ( 10 ) 

where Q3,Q4i,Q42 and Q43 are given in Appendix B. 
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2.3 Perturbation to the one-mode thermal state 



Assuming that a one-mode thermal state is perturbed, the perturbed characteristic function is 

x(p)=XT(jj)[l+s{a*li-aii*)], (11) 

where xt(p) — exp[— (N + i) |/i| 2 ] is the characteristic function of the one-mode thermal state, with N the average 

photon number of the state. The density operator can be obtained with p = J ^x(/i)D(— /i), where D(p) = exp^a* — 
p,*a] is the displacement operator, a^and a are the creation and annihilation operator of the system, respectively. By 
the method of integral within ordered product operator, the perturbed density operator is 

P = Pt + e(1 — v){aa) pt + a* pro)- (12) 

where pr = (1 — v) J^^Lo u ™ \ n ) ( n l * s ^ ne unperturbed one-mode thermal state with v = N / (N +1). It is clear that 
in the eigenbasis of pr the diagonal elements of the perturbation H = (1 — v){aa^ px + a* pto) is null, thus in the 
evaluation of Eq.® we only need to consider the off-diagonal elements of H. Evaluating Eq. ^ in the basis of pT,up 
to e 4 , for a null-diagonal H, the entropy will be 

S( P ) = S( PT ) - e 2 i g " m i log En + £3Qs ~ £ ^ Qil + Qi2 + Qi3 ^ + o(g4) (13) 



In our case, H nm = (1 — f)(ay / nE' n _i(5 n5m +i 4- a* E n ^fn + l<5 n;m _i). Thus Q3 = Q41 = due to the structure of H 7 
and Q43 = |a| 4 [ -^4-^ — h y^j log?;], Q42 = —2 |a| 4 [1 + jz^z log u]. The entropy of the perturbed state is 

S(p) = S(p T ) - 2e 2 M 2 log - - e 4 \af fi^- + \^ log -} + o(e 4 ). (14) 

3 The perturbation to the entropy of a degenerate system 
3.1 The entropy perturbation up to second order 

Note that Eq.© is an overall result regardless of the eigenvalue structure of the state po. The difference between the 
degenerate system and non-degenerate system comes when we evaluate the entropy. We now classify the eigenvectors 
according to the eigenvalues of po. Suppose the eigenvalue E n correspond to the eigenvector set {|n, rii) , \i = 1, . . .714 }, 
these eigenvectors span the subspace of dimension rid for E n . Define rid x trid matrix H nm with its entries being 

H nni<mmj = (n,ni\H\m,mj) . (15) 

Consider the first order derivative to the entropy, we have aa 

= E I" E "( E " + ty 2 TrH nn dt - }Z TrH nn \ogE n 

n ^ n 

= -J2TrH an logE n . (16) 

n 

Since TrH nn = Y^i=i ( n J n i \ H \ n > n i) is an invariant in the subspace, the first order perturbation to the entropy can 
be written as — J^k ^kk log-Efc, where k is the unified label for all the distinct eigenvectors (for some k' ^ k, we may 
have Ey = Thus the degenerate of the eigenvalues will not affect the expression of the first order perturbation of 
the entropy. With the notation of H nm , it follows the second order derivative to the entropy 

d 2 S(p ) _ v Tr[H m HU v Tr[H nm H^ m ] 

cfe 2 ^ E n 2- < E n -E m 8 " [ 1 

where Tr[-ff nn 7/J n ] = Tr[ff nn .ff* n ] = Ylij \Hnni,nn^\ is the summation of the absolute square of all entries of the 
matrix if nn ,not the summation of the absolute square of diagonal entries of the matrix H nn . 
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3.2 Perturbation to the two-mode thermal state 

Let pT — (1— v) Y^=o 1,71 \ n ) ( n \ ^ e single-mode thermal state, the direct product of pr with itself will result a two-mode 
P2 = Pt x Pt- The characteristic function of p 2 is x^ipi,^) = exp[— (N + + \p 2 \ )]■ We consider the perturbed 

characteristic function x(M1j M2) — X2(Mii A*2)[l + e(«Mi/ i 2 + a!*/-**/-*!)]- The perturbed density operator is p = P2 + sH, 
with 

H = (1 - v) 2 {a* a\a\p 2 + ap 2 aia 2 ), (18) 

where a\ and are the creation and annihilation operator for the two modes, respectively. The state P2 = (1 — 
v ) 2 J2^=o Sj=o w ™ l.?> n ~ ■?) 0') 71 — j\ ■ F° r any given n > 0, the state is (n -I- l)-fold degenerate. In the eigenbasis of p 2 , 
the entries of matrix H nm are 



(j,n- j\H \k,m- k) = (1 - v) 2 [a*8 jtk+1 5 n:m+2 E n - 2 Vj( n - j) 

+a6 jt k-iS n , m - 2 E n y/(J + l)(n -j + 1)], (19) 

where _E„ = (1 — w) 2 !)" is the eigenvalue of p 2 . Clearly, H nn now is a zero matrix, so that TrH nn = Tr[H an H^ n ] = 0. 
We only need to evaluate the second term of Eq. (jTTJ) . We have 

, E n -E m ^ = W 

" ^(i+l)(n-j + l)]logt;" 



2|a| 2 , 1 



log-. (20) 



Thus the entropy of the perturbed state will be 



1 + v v 



e 2 \a\\ 1 



S(p) = 2S(p T ) - *T"TT" l0 § - + °( £2 )- ( 21 ) 
1 + v u 

4 Density operator expansion 

In the calculation of the coherent information, even if we assume that the input state undergo a simple form perturbation 
p = po + eH, the output state £{p) and the combined output state (I ® £ ) |$) (^| may have a more complicated 
perturbation structures by the application of channel map £ and the purification from p to |\I>). So that we need to 
consider a more generic form of perturbation: 

p = po + + e 2 ^ (2) + e 3 H^ + ■■■ . (22) 

It is reasonable to require that TrH^ = is true for each n. The entropy perturbation can be obtained in due course, 
which is 

00 00 n „oq 

S(p) = Sipo)-J2 £nTr l H(n)l °SPo}-J2 en 12(- i y E / tiPo + ty'H^ipo + t)- 1 

xH^(po + ty 1 ■ ■ ■ H^ip + t)- 1 8 ll+l2+ ... +l ^ n . (23) 
We may express the second and the third order derivatives of the entropy in more explicit forms: 

k^^ 1 = - Tr r tipo + tr'H^ipo + ty'H^ipo + tyidt 

-Tr[H&> log po], (24) 
= Tr rtipo+ty^ipo+tr^dt 



-2Tr / tipo+t^H^ipo+ty^^ipo+t^dt 
Jo 

-Tr[H& log p ]. (25) 
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To verify the entropy expansion formula, we consider the displacement of thermal state, pd — D(ea)pTD^ (set). The 
characteristic function of the state is Xd{p) — Tr[pdD(pj\ = Tr[prD^ (ea)D(p,)D(ea)] — xt(p) exp[e(a*^ — oip*)]. 
Since the displacement operator is unitary, the entropy will not be changed by the displacement operation. Thus 
S{pd) = S(pr), the entropy S(pd) is not a function of e. If we expand the entropy of pd with respect to e, then the 
perturbation up to any order should be zero. Expanding with respect to e, the characteristic function is 

Xd {n) = Xt(m)[1 + £{a*p - ap*) + ^e 2 (a*p - ap*f + ■■■}. (26) 

The state is p d = p T + eH^ + e 2 H {2 ^ H , with 

ff (2) = —(1 — vf[a 2 a) 2 p T + PT a* 2 a 2 + 2 \a\ 2 a) p T a] - (1 - v) \a\ 2 p T . (27) 

The term has been treated in the former section, which give rise to the first term of the right hand side of Eq. 
(|24|) the value 2 \a\ 2 \ogv . And by direct calculation we have 

- Tr[H {2) log p ] = -2 \a\ 2 log v. (28) 
So d fgf ' 1 = 0. The first and the third derivatives of the entropy should be considering the symmetry of the system. 

5 Discussion and Conclusion 

In the examples presented, we obtain the density operator perturbation via the expansion of the characteristic function 
for continuous variable systems. This has some merits. First of all, since the characteristic function is a c-numbcr 
function, its expansion is simply Taylor expansion. We do not consider the Laurent expansion of the complex func- 
tion since for any state we have x(0) = 1. The characteristic function is defined as x(/i) = Tr[pD(p)], thus x(0) = 
Tr[pD(0)] = Tr[p] — 1. Secondly, we expand the characteristic function as x(a*) = Xo(j")[l £n S™=o c i,n-ifJ> 1 V*^] 

for the single-mode system (the multi-mode expansion can be obtained in due course). The unperturbed character- 
istic function Xo(/-0 is reserved as an integral nuclear for each order of perturbations. Thus the n— th perturbation 
= J fLi£^ (^) ^" =0 Ci^n-ip. 1 p*i D(—p) can be obtained with the difficulty of non-integrable. This form of expan- 
sion also implies that Tr[H^] = 0, which is not apparent at the first sight when we write down, for example, in 
Eq. ([27]) . although it can be proved by simple calculation. 

One of the obstacle that may be encountered in the entropy perturbation calculation is the eigenvalues of the 
unperturbed system. For instance, we consider the perturbation to a pure state, the unperturbed eigenvalues are 
and 1. This can be partially overcome by decomposing the perturbation matrix H as diagonal and off-diagonal part 
according to the eigensystem of the unperturbed density matrix po- Suppose p = po + s(Hq + H'), with Hq being 
diagonal and H' being off-diagonal in the eigenbasis of po- So po and Hq can be diagonalized simultaneously. Thus the 
perturbed state can be written as p = p' + eH' with a new expansion base state p' = po + eHq. Usually the eigenvalues 
of Po will differ from 0. 

We have derived the entropy perturbation formula for degenerate and non-degenerate system. Up to any give order 
of perturbation of the density operator, the entropy perturbation was written in a concise form via operator integral, 
the integral then can be evaluated in the eigenbasis of the unperturbed density operator of the system. 

Appendix A: The n-th order derivative of entropy 

^S^ 1 ==P 3 + P 4 , with 
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P 3 = -fm^f2 Tr PMP0t + ls)( l )(-ir- 1 

1=0 ^ ' 



p + le + t 



= h ^ b- 1 )' (1W rf f; E(-T) <i ¥- ( -T )< "¥ ( -T )m " B " Ey - o<,( - i)B 

2=0 ^ ' ^° m=n »i,—i n 

1=0 v 1 ' Ja 

POO 

= {-l) n n\ Trp Q (p +t)- 1 [H(p Q +t)- 1 } n dt. 
Jo 

The last equation cames from the identity E"=o( — -0' 0'" = (~l)" n - I n the fourth equation we have used the fact 
that Er=o( _1 )'0^ m = for a11 m < n - We now P rovc this tw0 results - Assuming = (1 - a;)" = ELot -1 )' 0^' 
let / ro (x) = = ELoM^O^' then /m(a) = Er=i4 m) ^ i ^ ) where 4 m) is the coefficient which 

obey the recursive relation c k m+V> = kc k m ^ + c k ™\. We always have cffl — 1. Note that = 0, for all k < n , and 

^ - (-l)»n!, for fc = n. Thus we have ££=<,(- l)'(i)i n = /"(I) = (-!)"«! and EzloM)'^" 1 = /-W = ° fOT 
all m < n. 

P4 can be obtained in the very similar way. 

Appendix B: The third and fourth order off-diagonal perturbation 

When the perturbation is off-diagonal in the eigenbasis of the unperturbed state, we have H nn — 0. The contribution 
to the third order perturbation should be 

poo 

Q 3 = / Tr{t( Po + ty 1 [H(p a +ty 1 f}dt. 
Jo 

Evaluating in the basis of p 0} we have 



at. 



n.m.k 



l (E n +t) 2 (E m +t)(E k +t) 



The summation should be taken for the case that n,m and k are all different since H nn = 0. After integrating and 
considering the permutation symmetry we have 



Q 3 = £ 



2Re(H nm H mk H kn ) Aog(E n /E k ) \og(E n /E m ). 



E m — E k (E n - E k ) E n — E„ 

n<m<k 

The contribution to the fourth order perturbation is 

/ T~f — T-f 1 T~fi 1, T~fi„~ 

-dt. 



nhj° (En+tnE m +t)(E l+ t)(E k +t) 1 



When all n, m, I, k are different with each other, the integral will be 



Qn — £ 2Re(H nm H m iHikH k n + H nTn H m iHikH k n + H nm H m iHikH kn ) 



n<m<l<k 



\og(E m /E n ) log{Ei/E n ) 



L {E m - E n )(E m - Ei)(E. m - E k ) (Ei - E n )(Ei - E m )(Ei - E k ) 
\og(E k /E n ) j 
(Ek — E n )(E k — E m )(E k — Ei) 
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where permutation symmetry are considered. When two of n,m,l,k are equal, that is (i) n = I 7^ m 7^ k and (ii) 
m = k 7^ n 7^ I, the integral will be 

p. \ f°° t \H nm H nk \ , \ f°° t \H mn H m i\ 

„tV° (E n + tr(E m + t)(E k +t) + n 2 ^ n J (E n + tf(E m + tf{E l+ t) ■ 

which is 

Q42 = (\HrunH n k\ 2 + \H mn H mk \ 2 + \H kn H km \ 2 ) 



n<m<k 



1 log(E m /E n ) \og{E k /E n ) 



E n (E n — E m )(E n — Ek) (E rn - E n ) 2 (E m — E k ) (E k — E n ) 2 (E k - E m ) 
When n — I and m = k, we have the integral 

t\H„ 



POO 



OO j. I TT 1 4 



+ t)3(25 rn +t)2 <ft ' 



which is 

v i rr i4r 1 / 1 1 \ 1 E n -i 

^43=2^1^1 i 2{E _ E \2^E~ + E~'~(E -E V ° g E~ V 

n<m ^yJ^n 1 -"m) J^n 1 - l m \*-<n Lj mj 
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